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ABSTRACT 



Single-star stellar population (ssSSP) models are usually used for spectral 

' r-|'' stellar population studies. However, more than 50% of stars are in binaries and 

evolve differently from single stars. This suggests that the effects of binary inter- 



P \ actions should be considered when modeling the stellar populations of galaxies 



c/3 \ and star clusters. Via a rapid spectral stellar population synthesis (RPS) model, 

we give detailed studies of the effects of binary interactions on the Lick indices 



and colours of stellar populations, and on the determination of the stellar ages 
^ I and metallicities of populations. 

\Q I Our results show that binary interactions make stellar populations less lumi- 

^ I nous, bluer, with larger age-sensitive Lick index (H/9) and smaller metallicity- 

sensitive indices (e.g., Mgb, Fe5270 and Fe5335) compared to ssSSPs. It also 
shows that when ssSSP models are used to determine the ages and metallicities 
Q I of stellar populations, smaller ages or metallicities will be obtained, when using 

two line indices (H/5 and [MgFe]) and two colours (e.g., u — R and R — K), 
respectively. Some relations for linking the stellar-population parameters ob- 



H I tained by ssSSPs to those obtained by binary-star stellar populations (bsSSPs) 

are presented in the work. This can help us to get some absolute values for 
stellar-population parameters and is useful for absolute studies. However, it is 
found that the relative luminosity-weighted stellar ages and metallicities obtained 
via ssSSPs and bsSSPs are similar. This suggests that ssSSPs can be used for 
most spectral stellar population studies, except in some special cases. 

Subject headings: galaxies: stellar content — galaxies: elliptical and lenticular, 
cD. 
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Introduction 



Stellar population synthe sis is a powerful technique to study the stellar contents of galax- 

Yungelson fc Tutukovlll997l . IXout et alJEoQTl . IPols et al ' 



ies and star cluster s (see, e.g 



1998 



Hurley et al.l 120071 ). It is also an important method to study the formation and evolution of 
galaxies. Simple stellar population (SSP) models that do not take binary interactions into ac- 
count ar e usually used for spectral stellar populat i on studies, as most models are ssSSP mod- 



els (e .g., iBruzual fc Charlotll2003l . IVazdekis 



19991. iFioc fcRocca-Volmerangd 119971. IWorthey 



1994f). However, as pointed by, e.g., iDuquennoy fc Mayod (ll99ll ). |Pinfield et al.l (120031 ). and 
Lodieu et al.l (120071 ). about 50% of stars are in binaries and they evolve differently from single 
stars. We can see this when comparing the isochrone of an ssSSP to that of a bsSSP (Fig.l). 
In fact, bsSSPs b etter fit the colour-magnitude diagrams (CMDs) of star clusters than ssSSPs 
( iLi fc Hanll2007l ). This suggests that binary interactions can affect stellar population synthe- 
sis studies significantly and it is, therefore, necessary to consider binary interactions. This 
i s also supported by some observational results, e.g., the Fa r-UV excess of ellipt i cal galaxies 



(iHan et al.l l2007f) and blue stragglers in star clusters (e.g., iDavies et al. 112004 iTian et al. 



20061 : kin et al.ll2007l ). These phenomena can be naturally explained via stellar populations 
with binary interactions, without any special assumptions. 

A few works have tried to model populations via binary stars and have presented some 
results o n the effects of binary in teractions on spectral stellar population synthesis. For 
example, IZhang et al.l (120041 . |2005| ) showed that binary interactions can make bsSSPs bluer 
than ssSSPs. However, there is not a more detailed investigation about how binary inter- 
actions affect the Lick indices and colours of stellar populatior is, and on t he de termination 
of stellar ages and metallicities. One of our previous works, iLi &: Hani (|2007l ). compared 
bsSSPs to ssSSPs, but various stellar population models are used. This makes it difficult to 
understand the effects of the changes of Lick indices and colours of populations that result 
only from binary interactions. Furthermore, it did not show how binary interactions affect 
the determinations of stellar ages and metallicities. In this case, we have no clear picture 
for the differences between the predictions of bsSSPs and ssSSPs, and do not well know the 
differences between luminosity-weighted stellar-population parameters (age and metallicity) 
determined by bsSSPs and ssSSPs. Because all galaxies contain some binaries, detailed 
studies of the effects of binary interactions on the Lick indices and colours of populations 
are important, as are the determinations of stellar-population parameters. In this work, we 
perform a detailed study of the effects of binary interactions on spectral stellar population 
synthe sis studies, via the rapid spectral population synthesis {RPS) model of iLi fc Han 

(EooJ- 



The paper is organized as follows. In Sect. 2, we briefly introduce the RPS model. 
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In Sect. 3, we study the effects of binary interactions on tlie isoclirones, spectral energy 
distributions (SEDs), Lick Observatory Image Dissector Scanner absorption line indices (Lick 
indices), and colours of stellar populations. In Sect. 4, we investigate the differences between 
stellar ages and metallicities fitted by ssSSPs and bsSSPs. Finally, in Sect. 5, we give our 
discussions and conclusions. 



2. The rapid spectral population synthesis model 



We take the results of the RPS model of iLi fc Had (120071 ) for the work, as there is 
no other available model. The RPS mod el calculated the evo lution of binaries and single 
stars via the rapid stellar evol ution code of iHurlev et al.l (120021) fhereaft e r Hur ley code) and 
took the spectral libraries of iMartins et al.l (120051 ) and IWestera et al.l (120021 ) (BaSeL 3.1) 
for spectral synthesis. The model calculated the high-resolution (0.3 A) SEDs, Lic k indices 
and co lour s for both bsSSPs and ssSSPs with the initial mass functions (IMFs) of ISalpeter 
(119551 ) and lChabrierl (120031 ). Note that the RPS m odel used a statistical isochrone database 
for modeling stellar populations (JLi &: Hani 120071 ). Each bsSSP contains about 50% stars 
that are in binaries with orbital periods less than 100 yr (the typical value of the Galaxy), 
and binary interactions such as mass transfer, mass accretion, common-envelope evolution, 
collisions, supernova kicks, angular momentum loss mechanism, and tidal interactions are 
considered when evolving binaries via Hurley code. Thus the RPS model is suitable for 
studying the effects of binary interactions on stellar population synthesis studies. However, 
some parameters such as the ones used for describing the common envelope prescription, 
mass-loss rates, and supernova kicks are free parameters and the default values in Hurley 
code, i.e., 0.5, 1.5, 1.0, 0.0, 0.001, 3.0, 190.0, 0.5, and 0.5, are taken in this work for wind 
velocity factor (/5w), Bondi-Hoyle wind accretion faction (ow), wind accretion efficiency factor 
(yUw), binary enhanced mass loss parameter (B^), fraction of accreted material retained in 
supernova eruption (e), common-envelope efficiency (oce), dispersion in the Maxwellian 
distribution for the supernovae kick speed (c3"k), Reimers coefficient for mass loss (77), and 
binding energy factor (A), respectively. These default values are taken because they have 
been tes ted by the develope r of Hurley code and seem more reliable. One can refer to the 
paper of iHurley et al.l (120021 ) for more details. In fact, many of these free parameters remain 
uncertain, and their uncertainties can possibly have great effect on our results. When we test 
the uncertainties caused by various oce and A, we find that the number of blue stragglers can 
be changed as large as 40% compared to the default case. However, it is extremely difficult to 
give detailed uncertainties in spectral stellar populatioi i synthesis due to t he free parameters, 
as we lack of constrains on these free parameters (see iHurley et al.l 120021 ) . Therefore, when 
we estimate the synthetic uncertainties in our RPS model, the uncertainties due to variation 
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of free parameters will not be taken into account. Because the f itted formulae used by Hurley 
code to evolve stars lead to uncertainties less than about 5% iHurley et al.l (120021 ) . we take 
5% as the uncertainties in the evolution of stars in the whole paper. The correctness of 
the results of the RPS model depends on how correct the default parameters of Hurley 
code are. In addition, the uncertainties in final generated spectrum caused by the spectral 
library and the method used for spectral synthesis are about 3% and 0.81%, respectively. 
Because a Monte Carlo technique is used by the RPS model to generate the star sample 
(2 000 000 binaries or 4 000 000 single stars in our work and it is twice of that of the model of 
Zhang et al.l 120041 ) of stellar populations, the number of stars can result in statistical errors 
in the Lick indices and colours of populati ons. According to our test, 1000 000 binaries is 
enough to get reliable Lick indices (see also lZhang et al.ll2005l ). but the near- infrared colours 
such as {I — K), {R — K), and (r — K) of old populations are affected by the Monte Carlo 
method. However, the errors caused by the Monte Carlo method are small, about 2% for a 
sample of 4000 000 stars. Note that in this work, an uniform distribution is used to generate 



the ratio (q, 0-1) of the mass of the secondary to that of the primary (JMazeh et al.l Il992 



Goldberg fc Mazehlll994l ). and then the mass of the secondary is calculated from that of the 
primary and q. The separation (a) of two components of a binary is generated following 
the assumption that the fraction of binary in an interval of log(a) is constant when a is 
big (lO-Rf;^ < a < 5.75 x IO^Rq) and it falls off smoothly when when a is small (< IORq) 



( iHan et al. Ill995l ). The distribution of a is written as 



a .p{a) 



^sep 



[a/aoY 



^sep; 



a < aQ 

ao < a < ai 



where Osep ~ 



0.070, ao = lOi^0, ai = 5.75 x IO^Rq and ip ^ 1.2. The eccentricity (e) of each 
binary s ystem using a unif orm distribution, in the range of 0-1, and e affects the results 
slightly JHurlev et aDl2002h . 



In addit i on, R PS model uses some methods different from those used in the work of 
Zhang et al.l (120041 ) to calculate the SEDs, Lick indices, and colours of populations. Be- 
sides RPS used a statis tical isochrone data base, it calculated the Lick indices directly from 
SEDs while the work of IZhang et al.l (120041 ) used some fitting formulae to compute the same 
in dices. The RPS ra odel calculated the colours of populations from SEDs, but the work 
of IZhang et al.l (2004r) co mputed colours by interpolating the photometry library of BaSeL 
1998 ). Furthermore, the RPS model used the more advanced version 



Leieune et al. 



2.2 



3.1 (IWestera et al. 



2OO2I ) of BaSeL library rather than version 2.2 to give the colours of the 



populations. The BaSeL 3.1 library overcomes the weakness of the BaSeL 2.2 library at low 
metallicity, because it has been colour-calibrated independently at all levels of metallicity. 
This makes the predictions of our model are more reliable. Another important point is that 
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the model of Zhang et al. did not present the near-infrared colours of stellar populations, but 
such colours are very important for disentangling the well-known age-metallicity degeneracy. 



3. Effects of binary interactions on stellar population synthesis 

3.1. The effects on isochrones of stellar populations 

The direct effect of binary interactions on stellar population synthesis is to change the 
isochrones of stellar populations, e.g., the distribution of stars in the surface gravity [log(f7)] 
versus effective temperature (Teg) grid, hereafter gT-grid. We investigate the di fferences 



between the isochrones of bsSSPs and ssSSPs. Stellar populations with the IMF of lSalpeter 



( 119551 ) are taken as our standard models for this work. The Salpeter IMF is actually not the 
best one for stellar population studies, although it is widely used. The reason is that the IMF 
is not valid for low masses. However, this IMF is reliable for stellar population synthesis, 
because low-mass stars contribute much less to the light of populations compared to high- 
mass stars. Some f urther studies will be given by taking some more reliable IMFs, e.g.. 



Kroupa et al.l (119931 ). Because the isochrone database used by this work divides the gT-grid 
into 1089 701 sub-grids, with 0.01 and 40K the intervals of log{g) and T^s, it is possible to 
compare the isochrones of bsSSPs and ssSSPs. The differences between the isochrones of two 
kinds of populations are calculated by subtracting the fraction of stars of bsSSP from that 
of its corresponding ssSSP, sub-grid by sub-grid. The ssSSP and its corresponding bsSSP 
have the same star sample, metallicity and age, and all their integrated specialties (SEDs, 
colours and Lick indices) are calculated via the same method. Therefore, the differences 
between the isochrones of a bsSSP and its corresponding ssSSP only result from binary 
interactions. For convenience, we call the difference a "discrepancy isochrone". Here we 
show the discrepancy isochrones for a few stellar populations in Figs. 2 and 3, for metal- 
poor {Z = 0.004) and solar-metalhcity {Z = 0.02) populations, respectively. Because it is 
found that the discrepancy isochrones of metal-rich [Z = 0.03) populations are similar to 
those of solar-metallicity populations, we do not show the results for metal-rich populations. 
Note that the results for populations with metallicities poorer than 0.004 are also given by 
our work, but we do not show them as the example for metal-poor populations as the RPS 
model did not give the SEDs and Lick indices for populations with metallicities poorer than 
0.004. This is actually limited by the spectral library used by the RPS model, which only 
supplies spectra for stars more metal rich than 0.002. As we see, some special stars, e.g., 
blue stragglers, are generated by binary interactions (see also Fig. 1). We also show that 
the differences between isochrones of old bsSSPs and ssSSPs are smaller than those of young 
populations, because the isochrones of old populations are dominated by low-mass stars, in 
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which binary interactions are much weaker. 



3.2. The effects on integrated features of populations 

The widely used integrated features of stellar populations are SEDs, Lick indices and 
colours. They are usually used for stellar population studies and are important. We inves- 
tigate the effects of binary interactions on them in this section. 



3.2.1. Spectral energy distributions 

To investigate how binary interactions affect the SEDs of stellar populations, we compare 
the SEDs of a bsSSP and an ssSSP that have the same age and metallicity. The difference 
between SEDs are simply called discrepancy SEDs. The absolute discrepancy SED for a 
pair of bsSSP and ssSSP is derived by subtracting the flux of the ssSSP from that of the 
bsSSP, as a function of wavelength. The discrepancy SEDs are mainly caused by blue 
stragglers and hot sub-dwarfs, as such stars are very hot and luminous. The changes of 
surface abundances of stars caused by binary interactions can also contribute to discrepancy 
SEDs. The absolute discrepancy SEDs for metal-poor {Z = 0.004), and solar-metallicity {Z = 
0.02) stellar populations are shown in Figs. 4, and 5, respectively. The absolute discrepancy 
SEDs such as those shown in the two figures can be easily used to add binary interactions into 
ssSSP models, but fractional discrepancy SEDs are more useful for understanding the effects 
of binary interactions. We show the fractional discrepancy SEDs of a few solar-metallicity 
populations in Fig. 6. 

As we see, binary interactions make stellar populations less luminous, but the flux 
in shortwave bands are changed by binary interactions weakly compared to that in long- 
wave bands. This mainly results from special stars generated by binary interactions, which 
contribute differently to flux in different bands. The differences between the SEDs of a 
bsSSP and its corresponding ssSSP decrease with increasing age or decreasing metallicity. 
In addition, it suggests that binary interactions can affect most Lick indices and colours of 
populations, because the flux changes caused by binary interactions are not 0, in the bands 
where widely used Lick indices and magnitudes are deflned. In this case, bsSSP and ssSSP 
models usually give different results for stellar population studies. Furthermore, because the 
effects of binary interactions on the SED flux of populations are about 11% on average, they 
are detectable for observations with spectral signal to noise ratio (SNR) grater than 10. In 
other words, the effects can be detected by most observations, as most reliable observations 
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Fig. 1. — Comparison of the shapes of the isochrones of a pair of solar-metalhcity {Z = 0.02) 
bsSSP and ssSSP. The figure is plotted by putting the isochrones of the bsSSP and ssSSP 
together. Black points show the isochrone of the ssSSP, and gray points the isochrone of the 
bsSSP. 
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Fig. 2. — Differences between the number distribution of stars in the isochrones of metal-poor 
[Z = 0.004) bsSSPs and ssSSPs.The darker the colour, the bigger the difference between star 
numbers of the bsSSP and ssSSP, (Nb - Ng). For each sub-grid, Nb is the number of stars in 
a bsSSP, while Ng the number of stars in an ssSSP and locating in the sub-grid. Note that 
the grid of log{g) versus T^s with log{g) range of -1.5 - 6 and T^fi range of 2 000 - 60 000 
K is divided into 1089 701 sub-grids, with intervals of 0.01 and 40 K for log{g) and Tes, 
respectively. A point in the figures corresponds to a sub-grid. 
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Fig. 3. — Similar to Fig.2, but for solar-metallicity {Z = 0.02) stellar populations. 
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Fig. 4. — Absolute differences between the SEDs of metal-poor {Z = 0.004) bsSSPs and 
ssSSPs, i.e., absolute discrepancy SEDs. The discrepancy SED of a pair of bsSSP and ssSSP 
is calculated by subtracting the flux of the ssSSP, Fg, from that of the bsSSP, Fb, as a 
function of wavelength. Each population contains 2 000 000 binaries or 4 000 000 single stars, 
whose initial mass follows the Salpeter IMF and with the lower and upper limits of 0.1 and 
100 Mq. 
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Fig. 5. — Similar to Fig. 4, but for solar-metallicity {Z = 0.02) stellar populations. 



-12- 



O - 




O 
I 

CM 

• 

O 
I 

ro 

o 

I 



Z=0.02, *=1.5 Gyr 



;i ^^-^-''^^^^^^--v^^ww-^.-,^ 



Z=0.02, t=5 Gyr 




'- (d) 



V— *— «iiii<i» 



Z=0.02, t=12 Gyr 



1 I I I I L. 



± 



3000 



4000 



5000 
X(A) 



6000 



7000 



Fig. 6. — Fractional discrepancy SEDs of solar-metallicity stellar populations. Symbols Fg 
and Fb have the same meanings as in Fig. 4. 
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have SNRs greater than 10. 



3.2.2. Lick indices 



Lick indices are the most widely used indices in stellar population studi es, because the y 
can disentangle the well-known stellar age-metallicity degeneracy (see, e.g.. IWortheylll994l ). 
If binary interactions are taken into account, some results different from those determined 
via ssSSP models will be obtained, which was suggested by the study of the differences 
between the SEDs of bsSSPs and ssSSPs. Here we test how binary interactions change the 
Lick indices of stellar populations when comparing to those of ssSSPs. In Fig. 7, we show 
the differences between four widely used indices of bsSSPs and those of ssSSPs. 



The indices are calculated from SEDs on the Lick system (jWorthey et al.lll994l ) directly. 
As we see. Fig. 7 shows that binary interactions make the H/? index of a population larger 
by about 0.15 A, while making Mgb index smaller by about 0.06 A and Fe indices smaller 
by more than about 0.1 A, compared to ssSSPs. Therefore, the changes in Lick indices are 
usually larger than typical observational uncertai nties (about 0.07 A f or H/5 index and 0.04 
A for metal-line indices according to the data of iThomas et al.ll2005l ). For fixed metallic- 
ity, the effects of binary interactions on both age- and metallicity-sensitive indices become 
stronger with increasing age when stellar age is less than about 1.5 ~ 2 Gyr, and then the 
effects become weaker with increasing age. The reason is that binary interactions change the 
isochrones of populations most strongly near 1.5 ~ 2 Gyr as the first mass transfer between 
two components of binaries peak and a lot of blue stragglers are generated near 1.5 ~ 2 Gyr 
according to the star sample of bsSSPs, and the light of old populations is dominated by 
low-mass binaries. The interactions between two components of low-mass binaries are usu- 
ally weaker than for high-mass binaries. The effects of binary interactions on the isochrones 
are tested quantitatively using the numbers of stars with log(5') < 4.0 and log(Teff) > 3.75, 
because these stars are very luminous and contribute a lot to the light of their populations. 
Our result shows that binary interactions change the number distribution of stars in the 
above log{g) and log(Teff) ranges most significantly when stellar age is from 1.5 to 2. In 
addition, from Fig. 7, we find that for a fixed age, binary interactions affect the H/3 and 
Fe5270 indices of metal-poor populations more strongly than those of metal-rich populations, 
while they affect the Mgb and Fe5335 indices of metal-poor populations more weakly. As 
a whole, using ssSSPs and bsSSPs, various ages and metallicities will be measured for the 
same stellar population, via popular Lick-index methods such as the H/5 & [MgFe] method, 
which determines the ages and metalli cities of populations by comparing the observational 
and theoretical H/5 and [MgFe] indices (IThomas et al.ll2003l ). Note that the evolution of the 
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Fig. 7. — Differences between four widely used Lick indices of bsSSPs and ssSSPs. Tfie 
difference in a Lick index is calculated by subtracting tlie value of a ssSSP from that of its 
corresponding bsSSP, wliich has the same age and metalhcity as the ssSSP. Star, circle, cross, 
and square are for populations with metalhcities of 0.004, 0.01, 0.02, and 0.03, respectively. 
Note that the differences are averaged in each bin. 
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differences between the Lick indices of ssSSPs and bsSSPs were not shown before. 

3.2.3. Colour indices 
Beca use colours are u seful for estimating the ages and metallicities of distant galaxies 



[see, e.g. 



Li fc Harul2008l ). we investigate the effects of binary interactions on them. We 
use a method similar to that used for studying the Lick indices of stellar populations. Our 
detailed results are shown in Fig. 8. In the figure, the differences between two UBVRIJHK 
colours, {B — V) and {B — K)., a ugriz colour on the photometric system used by Sloan 
Digital Sky Survey (hereafter SDSS system), (m — r), and a composite colour, {r — K), of 
bsSSPs and those of ssSSPs are shown, respectively. Note that the (r — K) colour consists 
of a Johnson system magnitude, K, and an SDSS system magnitude, r. As we see, for a 
fixed age and metallicity, binary interactions make the colours of most stellar populations 
bluer than those predicted by ssSSPs. This mainly results from blue stragglers generated by 
binary interactions, because such stars are very luminous and blue. It suggests that we will 
get different stellar metallicities and ages for galaxies via bsSSP and ssSSP models, using a 
photometric method. When comparing to the typical colour uncertainties (0.12, 0.06, 0.13, 
0.10, and 0.01 mag for B, V, K, u, and r magnitudes, respectively), the changes [e.g., about 
-0.04, -0.15, and -0.08 mag for (B — V), (B—K), and (u—r), respectively] of colours caused by 
binary interactions are similar to, but somewhat less than, typical observational errors. Note 
that the photometric uncertainties are estimated using the data of some publications, SDSS, 
and Two Micron All Sky Survey (2MASS). The uncertainties actually depend on surveys. 
The observational uncertainties of u and K magnitudes may be smaller when taking the data 
of other surveys instead of those of SDSS and 2MASS. In addition, similar to Lick indices, 
the differences between the colours of two kinds of stellar populations peak near 2 Gyr. 



4. The effects of binary interactions on the determination of 
stellar-population parameters 

Two stellar-population parameters, i.e., stellar age and metallicity, are crucial in the 
investigations of the formation and evolution of galaxies. We investigate the effects of binary 
interactions on the estimates of the two parameters. We try to fit the stellar-population 
parameters of bsSSPs with various ages and metallicities using ssSSPs, via both Lick-index 
and photometric methods. Because observations show that about 50% of stars in the Galaxy 
are in binaries, bsSSPs should be more similar to the real stellar populations of galaxies and 
star clusters. Therefore, the stellar-population parameters fitted via ssSSPs (hereafter ss- 
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Fig. 8. — Differences between four colours of bsSSPs and ssSSPs. The differences in a colour 
are calculated by subtracting the colour of a ssSSP from that of its corresponding bsSSP 
(with the same age and metallicity). Symbols have the same meanings as in Fig. 6. The 
colour {u — r) is on the SDSS-ugriz system and (r — K) is a composite colour that consists 
of a Johnson magnitude (K) and a SDSS-ugriz magnitude (r). The differences are averaged 
in each bin. 
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fitted results, represented by tg and Zg) should be different from the results obtained via 
bsSSPs (bs-fitted results, tb and Zb). The detailed differences are shown in this section. 



4.1. Lick-index method 

In a widely used method, i.e.. Lick-index method, we fit the stellar ages and metallicities 
of populations by two indices, i.e., H/3 and [MgFe] = v/Mgb.(0.72Fe5270 + 0.28Fe5335), after 



Thomas et al.l (120031 ). Thus the results are slightly affected by a-enhancement and stellar 
population models (e.g., our RPS model) without any a-enhancement compared to the sun 
can be used to measure stellar-population parameters. The differences between bs- and 
ss-fitted stellar-populations parameters of populations with four metallicities (0.004, 0.01, 
0.02, and 0.03) and 150 ages (from 0.1 to 15Gyr) are tested. In the test, we try to fit 
the stellar ages and metallicities of testing bsSSPs via an H/3 versus [MgFe] grid of ssSSPs. 
Because ssSSPs predict different Lick indices for populations compared to bsSSPs, when we 
use ssSSPs to fit the ages and metallicities of our testing bsSSPs, the results obtained are 
different from the real parameters of bsSSPs, i.e., bs-fitted parameters. From an H/5 versus 
[MgFe] grid of ssSSPs (Fig. 9), we can see this clearly. In detail, the ss-fitted stellar ages are 
less than the bs-fitted ones, by ~ 5 Gyr. The maximal difference is larger than the typical 
uncertainty (< 2 Gyr) in stellar population studies (see Fig. 9). The older the populations, 
the bigger the difference between ages fitted via bsSSPs and ssSSPs, although the differences 
between the Lick indices of old bsSSPs and ssSSPs are smaller (see Section 3). The reason 
is that the differences among the Lick indices of populations with different ages are much 
less for old populations than for young populations (see Fig. 9 for comparison). Therefore, 
the ss-fitted ages of galaxies can be much less than bs-fitted ages, because most galaxies, 
especially early-type o nes, have old (7~8 G yr) populations and their metallicities are not 



big (peak near 0.002) (iGallazzi et al.ll2005l ). However, ss-fitted metallicities of populations 
are similar to bs-fitted values, compared to the typical uncertainties (~ 0.002). Therefore, 
if some stars are binaries, smaller ages will be measured via comparing the observational 
H/3 and [MgFe] indices of galaxies with those of theoretical ssSSPs. This is more significant 
for metal-poor stellar populations. In our testing populations, on average, the ss-fitted 
metalhcities are 0.0010 poorer than bs-fitted values, while ssSSP fitted ages are younger 
than bs-fitted values, by 0.3 Gyr for all and 1.8 Gyr for old (> 7 Gyr) testing populations. 
In this work, the ss-fitted stellar ages and metallicities of testing bsSSPs are obtained by 
finding the best-fit populations in a grid of theoretical populations with intervals of stellar 
age and metallicity of 0.1 Gyr and 0.0001, respectively. A least square method is used in the 
fit. In addition, it is found that the bs-fitted ages of populations can be calculated from the 
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20031 ). Dashed and solid lines are for constant metaUicity and constant age, respectively. Er 
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ss- fitted ages and metallicities (with RMS of 1.45 Gyr) via the equation 

tb = (0.17 + 8.27Z,) + (1.38 - UA5Z,%, (2) 

where tb, Zs and tg are the bs-fitted age, ss-fitted metalhcity and age, respectively. The 
relation between bs-fitted ages and ss-fitted stellar-population parameters can be seen in 
Fig. 10. We find that ss-fitted ages are usually less than the bs-fitted ages of populations, 
and the poorer the metallicity, the larger the differences between the bs- and ss-fitted ages. 
Note that Eq. (2) is not very accurate for metal-poor {Z = 0.004) and old (age > 11 Gyr) 
populations. The reason is that the B.j3 index increases with age for metal-poor and old 
populations, while it decreases with age for other populations. 

4.2. Photometric method 

In the photometric method, we fit stellar-population parameters respectively via two 
pairs of colours, i.e., [{u — R), {R — K)] and [{u — r), (r — K)]. The two pairs are shown 
to have the ability to constrain the ages and metallicities o f populations an d can be used to 



study the stellar populations of some distant galaxies (see lLi fc Hanll2008l ). The test shows 
that ss-fitted metallicities are poorer than the bs-fitted metallicities of populations. When 
taking [{u — R), {R — K)] for this work, on average, ss-fitted metallicities are 0.003 smaller 
than bs-fitted values. It is 0.0035 when taking the pair [{u — r), (r — K)]. The distribution 
of a few testing bsSSPs in the {u — R) versus {R — K) grid of ssSSPs is shown in Fig. 
11. In particular, it is found that the ss-fitted ages are correlated with the bs-fitted ages 
of populations, which is independent of metallicity. The relation (with a RMS of 0.72 Gyr) 
between ss- and bs-fitted ages of populations can be written as 

tb = 0.24 + 0.93ts, (3) 

where tb and t^ are bs- and ss-fitted ages, respectively. It shows that the bs- and ss-fitted 
ages are similar. The equation is clearly different from Eq. (2), because colours are usually 
less sensitive to metallicity compared with metal-line Lick indices. The relation between bs- 
and ss-fitted ages of populations is shown in Fig. 12. The figure shows the approximate 
relation between the bs- and ss-fitted ages of populations, which is nearly independent of 
metallicity. The equation is possibly useful to estimate the absolute ages of distant galaxies 
and star clusters. Note that the relation is presented for populations younge r than 14 Gyr, 



becau se the age of the universe is shown to be smaller than about 14 Gyr (JBennett et al. 



20031). 



When we use {u — r) and (r — K) colours to estimate the stellar-population parameters 
of populations, we find that ss-fitted metallicities are about 0.0035 smaller than the bs-fitted 
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Fig. 10. — The relation between bs- fitte d ages of bsSSPs and tlieir ss-fitted stellar-population 
parameters, when using H/3 and [MgFe] (JThomas et al.ll2003l ) to measure the ages and metal- 
licities of testing bsSSPs. Star, circle, cross, and square point represent stellar populations 
with real metallicities of 0.004, 0.01, 0.02, and 0.03, respectively. Sohd, dash-dotted, dotted, 
and dashed lines show the fitted relations between the bs-fitted ages and ss-fitted results of 
populations, for ss-fitted metallicities of 0.003, 0.009, 0.019, and 0.029, respectively. 
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Fig. 11. — Some bsSSPs on the {u — K) versus {R — K) grid of ssSSPs. The position "b" 
shows the bs-fitted age and metalhcity of a bsSSP, and "s" shows the ss-fitted values. Dashed 
and sohd hues are for constant metalhcity and constant age, respectively. Error bars give 
typical observational uncertainties taken from the NED, SDSS and 2MASS surveys. 
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Fig. 12. — The relation between bs- and ss-fitted ages and metallicities of populations, when 
using {u — R) and {R — K) colours to measure the parameters. Points are for populations 
with Salpeter IMF and have the same meanings as in Fig. 10. Sohd and dashed hues show 
the fitted relations between bs- and ss-fitted ages, for populations with Salpeter and Chabrier 
IMFs, respectively. 
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values. The bs- and ss-fitted ages of populations can be approximately transformed by 

tb = 0.28 + 0.91ts, (4) 

where tb and tg are the bs- and ss-fitted ages, respectively. The RMS of the fitted relation is 
1.00 Gyr. The equation is similar to Eq. (3) but with larger scatter. This results from the 
various metallicity and age sensitivities of the colours used. One can see Fig. 13 for more 
details about the relation. 

As a whole, from both the results obtained by Lick-index and photometric methods, we 
are shown that bs-fitted stellar-population parameters increase with ss-fitted ones. Therefore, 
using bsSSP models instead of ssSSP models, similar results for relative studies of stellar- 
population parameters of galaxies will be obtained. However, if one wants to get the absolute 
stellar-population parameters of galaxies and star clusters, the effects of binary interactions 
should be taken into account, especially for metal-poor populations. It can be conveniently 
done by taking the average metallicity deviations and the relations between the bs-fitted 
ages and ss-fitted results of populations, which were shown above. 

4.3. Results for populations with Chabrier initial mass function 



Some stellar populations with Salpeter IMF (ISalpeterl Il955l ) were taken as standard 
models for the work, but even if some populations with other IMFs were taken instead, 
we can obtain similar results. We have a test using populations with the Chabrier IMF 



(IChabrierll2003l ) . The result shows that ss-fitted metallicities are 0.0011 less on average than 
bs-fitted results, when taking H/? and [MgFe] for measuring stellar-population parameters. 
The bs-fitted ages and ss-fitted stellar-population parameters have a relation of tb = (—0.06 + 
20.63Zs) + (1.46— 18.76Zs)ts, where Zg is the ss-fitted metalhcity, while tb and 4 the bs- and ss- 
fitted ages, respectively. When we used {u—R) and (R—K) to estimate the stellar-population 
parameters of populations, ss-fitted metallicities were shown to be 0.0031 smaller than bs- 
fitted values, and bs-fitted ages can be calculated from ss-fitted results via tb = 0.46 + 0.89ts. 
A similar relation for the results fitted by [u — r) and (r — K) is tb = 0.41 + O.SStg, with a 
deviation of 0.0039 in metallicity. As a whole, the relations between bs- and ss-fitted results 
obtained via populations with Salpeter and Chabrier IMFs are similar, compared to the 
typical uncertainties of stellar-population parameter studies. The comparisons of the results 
obtained via two IMFs can be seen in Figs. 11 and 12. 
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5. Discussions and Conclusions 

We investigated the effects of binary interactions on tlie isoclirones, SEDs, Lick indices 
and colours of simple stellar populations, and then on the determination of light-weighted 
stellar ages and metallicities. The results showed that binary interactions can affect stellar 
population synthesis studies significantly. In detail, binary interactions make stellar popula- 
tions less luminous and bluer, while making the H/3 index larger and metal line indices smaller 
compared to ssSSPs. The colour changes (2~5%) caused by binary interactions are smaller 
than the systematic errors (about 6%) of the RPS model while similar to observational errors 
(4~7%). Note that the systematic error of 6% did not take the uncertainties due to the free 
parameters of the star model into account (see Sect. 2). The changes (3~6%) of Lick indices 
caused by binary interactions are somewhat smaller than the systematic errors (about 6%) of 
stellar population synthesis model but larger than observational errors (1~4%). Therefore, 
if we measure luminosity-weighted stellar-population parameters (metallicity and age) via 
bsSSPs instead of ssSSPs, higher (0.0010 on average) metallicities and significantly larger 
ages will be obtained via a Lick-index method, and significantly higher (about 0.0030) metal- 
licities and similar ages will be obtained via a photometric method. Because simple stellar 
population models are usually used for studying the populations of early-type galaxies or 
globular clusters, which possibly have old (> 7 Gyr) and relative metal-poor populations, the 
changes (~ 1.8 Gyr in age and 0.0030 in metallicity) caused by binary interactions in stellar 
ages and metallicities are larger than the typical uncertainties. In particular, we found that 
the relative results of stellar population studies obtained by ssSSPs and bsSSPs are similar. 
The bs-fitted stellar-population parameters can be calculated from the ss-fitted ones, via 
equations presented by the paper. 

The relations between bs-fitted ages and ss-fitted stellar-population parameters are use- 
ful for some special investigations. For example, when studying the age of the universe via 
the stellar ages of some distant globular clusters, we can estimate the absolute age of star 
clusters using ss-fitted results. Although the results shown in this paper can help us to give 
some estimates for the absolute stellar ages and metallicities of galaxies, it is far from getting 



accur ate values because of the large uncertainties in stellar population models (see, e.g., lYi 



20031 ). In addition, different stellar population models usually give different absolute results 



for stellar population studies. Note that the results obtained by the Lick-index m ethod are 



affected slightly by a-enhancement, according to the work of iThomas et al.l (120031 ) . but this 
is not the case for the results obtained by photometric methods. 

In this work, all bsSSPs contains about 50% binaries with orbital periods less than 100 
yr (the typical value of the Galaxy). If the binary (with orbital periods less than 100 yr) 
fraction of galaxies are different from 50%, the results shown in the paper will change. The 
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higher the fraction of binaries, the larger the difference between ss- and bs-fitted stellar- 
population parameters. Thus the results obtained by in this paper may not be proper for 
investigating galaxies or star clusters with binary fractions obviously different from 50%. 
Furthermore, when bui lding bsSSPs, w e assumed that the masses of the two components of 



a binary are correlated (JLi fc Hanll2007l ). according to previous works. We did not try to take 
other distributions for secondary mass and binary period in this work, because it is limited 
by our present computing ability. We will give further studies in the future. The differences 
between the Lick indices and colours of bsSSPs and ssSSPs do not evolve smoothly with 
age. This possibly relates to the method used to calculate the integrated features of stellar 
populations. In fact, the Monte Carlo method usually leads to some scatter (about 2%) in 
the integrated features of populations. The analytic fits and the binary algorithm used by 
Hurley code can also lead to some scatter (about 5%). 

We investigated the effects of binary interactions via only some simple stellar popula- 
tions, but the real populations of galaxies and star clusters are usually not so simple. In 
other words, the populations of galaxies and star clusters seem to be co mposite stellar pop- 



ulations including populations with various ages and metallicities (e.g., lYi et al.ll2005l ). It 
seems that the effects of binary interactions and population-mixing are degenerate. This is 
a complicated subject, which requires further study. 
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the Youth Foundation of Knowledge Innovation Project of The Chinese Academy of Sciences 
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